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DESCRIPTION 

ULTRASONIC APPARATUS AND METHOD FOR MEASURING THE 
CONCENTRATION AND FLOW RATE OF GAS 
5 Technical Field 

The invention relates to ultrasonic apparatus and 
method for measuring the concentration of oxygen gas in a 
sample gas and flow rate of the sample gas, which is 
supplied from an oxygen concentrator used for a medical 
10 purpose. 

Background Art 

It is well known that the propagation velocity of 

ultrasonic waves through a sample gas is presented by a 
15 function of the concentration and the temperature of the 

sample gas. The velocity C(m/sec) of ultrasonic waves 

propagating through a stationary gas is presented by 

flowing equation ( 1 ) with mean molecular weight M and the 

temperature T ( K ) . 
20 C=(KRT/M) 1/2 ...(1) 

Where: 

k: ratio of molecular specific heat at constant volume 

and molecular specific heat at constant pressure 
R: gas constant 

25 Therefore, measuring the velocity of ultrasonic 

waves C(m/sec) propagating through a sample gas and the 
temperature T(K) of the sample gas will provide the mean 
molecular weight M of the sample gas through a 
calculation. For example, the mean molecular weight M of 

30 a sample gas containing an oxygen-nitrogen gas mixture of 
a mixture ratio P: ( 1-P ) ( O^P^l ) will be calculated by 
following equation (2). 

M=M 02 P+M N2 (1-P) ...(2) 

Where : 

35 M Q2 : Molecular Weight of oxygen gas 

M N2 : Molecular Weight of nitrogen gas 
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Therefore, the oxygen concentration P will be 
obtained through a calculation on the basis of the 
measurement of mean molecular weight M. When the sample 
gas is an oxygen-nitrogen mixture, k=1.4 is reasonable 
5 over a wide range of the oxygen-nitrogen mixture ratio. 

When the velocity of ultrasonic waves propagating 
through a sample gas is C(m/sec) and the flow velocity of 
the sample gas is V(m/sec), the velocity of ultrasonic 
waves Cxfm/sec) propagating in the forward direction 
10 relative to the sample gas flow is C X =C+V, and the 

velocity of ultrasonic waves C 2 (m/sec) propagating in the 
backward direction relative to the sample gas flow is 
C 2 =C-V. Therefore , the velocity of the sample gas flow 
V(m/sec) is calculated by following equation (3). 
15 V=(C 1 -C 2 )/2 ...(3) 

The flow rate (mVsec) of the sample gas will be 
obtained by multiplying the velocity of the sample gas 
flow by the sectional area (m 2 ) of the conduit through 
which the sample gas flows. 
20 Methods and apparatuses for measuring the 

concentration of a certain gas or the flow velocity of a 
sample gas, by using the above principle, on the basis of 
the propagation velocity or the propagation time of 
ultrasonic waves through the sample gas have been 
25 developed. For example, Japanese Unexamined Patent 

Publication (Kokai) No. 6-213877 describes an apparatus 
for measuring the concentration and the flow rate of a 
sample gas by measuring the propagation time of 
ultrasonic waves propagating between two ultrasonic 
30 transducers opposingly disposed in a conduit through 

which the sample gas flows. Further, Japanese Unexamined 
Patent Publications (Kokai) No. 7-209265 and No. 8-233718 
describe an apparatus for measuring the concentration of 
a certain gas contained in a sample gas by measuring the 
35 propagation velocity or propagation time of ultrasonic 

waves propagating through a volume with a reflecting type 
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apparatus including an ultrasonic transducer and an 
opposingly disposed reflector. 

in such a method and an apparatus for measuring the 
concentration and the flow rate by using the propagation 
5 velocity of the ultrasonic waves, it is necessary to 

accurately measure the propagation time of the ultrasonic 
waves. However, the signal generated on the basis of the 
received ultrasonic waves always includes noise 
component, which makes difficult to determine the moment 

10 when ultrasonic waves are received by the ultrasonic 
transducer. Therefore, the propagation time of 
ultrasonic waves is indirectly estimated through a 
complex signal processing procedure or a complex 
hardware. For example, Japanese Unexamined Patent 

15 Publication (Kokai) No. 9-318644 describes a method for 

measuring a propagation time of ultrasonic waves in which 
the waveform of the received ultrasonic waves is 
integrated. After the results of the integration of the 
waveform reach a predetermined vale, the first zero-cross 

2 0 time instant is determined as the propagation time of the 
ultrasonic waves for the measurement of the flow rate. 
According to the method, the timing of the generation of 
the zero-cross signal is not fluctuated even if the 
amplitude of the received waves is fluctuated to some 

25 extent. Therefore, the obtained zero-cross time instant 
is relatively close to the moment when the ultrasonic 
waves actually reach. However, the obtained zero-cross 
time instant is not real propagation time of the 
ultrasonic waves. In particular, when the concentration 

30 is measured, the measurement error is strongly affected 
by the difference between the real propagation time and 
the zero-cross time instant. 

Further, Japanese Unexamined Patent Publication 
(Kokai) No. 60-138422 describes a flow rate measuring 

35 device in which an envelope curve is calculated on the 
basis of the waveform of the received ultrasonic waves. 
The rise time of the envelope curve is calculated by an 
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approximate equation to estimate the ultrasonic 
propagation time. However, a hardware is necessary to 
sample the received ultrasonic waves and a complex signal 
processing is necessary to calculate the envelope curve 
5 based on the sampled waveform. Therefore , according to 
the invention of JPP '422, it is difficult to provide a 
compact device with low cost. 

Disclosure of the invention 

10 The object of the invention is to provide an 

ultrasonic apparatus and method for measuring 
concentration and flow rate of gas, which allow to 
accurately measure the concentration and flow rate of a 
sample gas without a complex signal processing and an 

15 additional hardware. 

According to the present invention, there is 
provided an ultrasonic apparatus for measuring the 
concentration and flow rate of a sample gas, comprising: 
a conduit for flowing the sample gas; 

20 a first ultrasonic transmission-reception device 

mounted to the inside of the conduit; 

a second ultrasonic transmission-reception device 
mounted to the inside of the conduit to face the first 
ultrasonic transmission-reception device; 

25 a transmission-reception switch for switching the 

operation mode of the first and second ultrasonic 
transmission-reception devices between a transmission 
mode for transmitting ultrasonic waves and a reception 
mode for receiving ultrasonic waves; 

30 a temperature sensor, disposed in the conduit, for 

measuring the temperature of the sample gas flowing 
through the conduit; 

the first ultrasonic transmission-reception device 
generating forward ultrasonic waves relative to the flow 

35 direction of the sample gas when the device is in the 

transmission mode and generating backward waveform when 
the device is in the reception mode on the basis of the 
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received ultrasonic waves generated by the second 
ultrasonic transmission-reception device; 

the second ultrasonic transmission-reception device 
generating backward ultrasonic waves relative to the flow 
5 direction of the sample gas when the device is in the 
transmission mode and generating forward waveform when 
the device is in the reception mode on the basis of the 
received ultrasonic waves generated by the first 
ultrasonic transmission-reception device; 
10 means for generating trigger signals when the 

forward and backward waveforms pass over a predetermined 
level; 

means for generating forward and backward zero-cross 
signals when the forward and backward waveforms pass over 

15 a zero level; 

propagation time calculation means, coupled to the 
temperature sensor, the trigger signal generating means 
and the zero-cross signal generating means, for (1) 
calculating a possible propagation time range on the 

20 basis of the gas temperature detected by the temperature 
sensor, (2) determining whether or not the phases at 
which two first trigger signals, respectively generated 
on the basis of the forward and backward waveforms, 
coincide with each other, (3) processing the zero-cross 

25 signals so that the phases coincide with each other if 
they do not coincide with each other, (4) obtaining 
reference zero-cross time instant by calculating mean 
value of the forward and backward zero-cross time 
instants, (5) obtaining an ultrasonic reception point by 

30 subtracting an integral multiple of the cycle of the 

ultrasonic waves so that the results of the subtraction 
falls into the possible propagation time range and (6) 
estimating the ultrasonic propagation time on the basis 
of the ultrasonic reception point, 

35 Further, according to another feature of the 

invention, there is provided a method of measuring the 
concentration of sample gas flowing through a conduit, 
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comprising the steps of: 

generating forward ultrasonic waves relative to the 
flow direction of the sample gas; 

generating backward ultrasonic waves relative to the 
5 flow direction of the sample gas; 

measuring the temperature of the sample gas flowing 
through the conduit; 

generating trigger signals when the forward and 
backward waveforms pass over a predetermined level; 
10 generating forward and backward zero-cross signals 

when the forward and backward waveforms pass over a zero 
level; 

calculating a possible propagation time range on the 
basis of the gas temperature detected by the temperature 
15 sensor; 

determining whether or not the phases at which two 
first trigger signals, respectively generated on the 
basis of the forward and backward waveforms, coincide 
with each other; 

20 processing the zero-cross signals so that the phases 

coincide with each other if they do not coincide with 
each other; 

obtaining reference zero-cross time instant by 
calculating mean value of the forward and backward zero- 

25 cross time instants; 

obtaining an ultrasonic reception point by 
subtracting an integral multiple of the cycle of the 
ultrasonic waves so that the results of the subtraction 
falls into the possible propagation time range; and 

30 estimating the ultrasonic propagation time on the 

basis of the ultrasonic reception point. 

Further, according to another feature of the 
invention, there is provided an oxygen concentration 
system for generating an oxygen enriched gas, comprising 

35 an oxygen concentration apparatus for generating an 

oxygen enriched gas by adsorbing nitrogen to remove the 
nitrogen from the air; and 
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an ultrasonic apparatus for measuring the 
concentration of the oxygen in the oxygen enriched gas 
and flow rate of the oxygen enriched gas, the ultrasonic 
apparatus comprising: 
5 a conduit for receiving and flowing the oxygen 

enriched gas; 

a first ultrasonic transmission-reception device 
mounted to the inside of the conduit; 

a second ultrasonic transmission-reception device 
10 mounted to the inside of the conduit to face the first 
ultrasonic transmission-reception device; 

a transmission-reception switch for switching the 
operation mode of the first and second ultrasonic 
transmission-reception devices between a transmission 
15 mode for transmitting ultrasonic waves and a reception 
mode for receiving ultrasonic waves; 

a temperature sensor, disposed in the conduit, for 
measuring the temperature of the oxygen enriched gas 
flowing through the conduit; 
20 the first ultrasonic transmission-reception device 

generating forward ultrasonic waves relative to the flow 
direction of the oxygen enriched gas when the device is 
in the transmission mode and generating backward waveform 
when the device is in the reception mode on the basis of 
25 the received ultrasonic waves generated by the second 
ultrasonic transmission-reception device ; 

the second ultrasonic transmission-reception device 
generating backward ultrasonic waves relative to the flow 
direction of the oxygen enriched gas when the device is 
30 in the transmission mode and generating forward waveform 
when the device is in the reception mode on the basis of 
the received ultrasonic waves generated by the first 
ultrasonic transmission-reception device; 

means for generating trigger signals when the 
35 forward and backward waveforms pass over a predetermined 
level; 

means for generating forward and backward zero-cross 
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signals when the forward and backward waveforms pass over 
a zero level; 

propagation time calculation means, coupled to the 
temperature sensor, the trigger signal generating means 
5 and the zero-cross signal generating means, for (1) 
calculating a possible propagation time range on the 
basis of the gas temperature detected by the temperature 
sensor, (2) determining whether or not the phases at 
which two first trigger signals, respectively generated 

10 on the basis of the forward and backward waveforms, 

coincide with each other, (3) processing the zero-cross 
signals so that the phases coincide with each other if 
they do not coincide with each other, (4) obtaining 
reference zero-cross time instant by calculating mean 

15 value of the forward and backward zero-cross time 

instants, (5) obtaining an ultrasonic reception point by 
subtracting an integral multiple of the cycle of the 
ultrasonic waves so that the results of the subtraction 
falls into the possible propagation time range and (6) 

20 estimating the ultrasonic propagation time on the basis 
of the ultrasonic reception point. 

Further, according to another feature of the 
invention, there is provided an oxygen concentration 
system for generating an oxygen enriched gas, comprising: 

25 an oxygen concentration apparatus for generating an 

oxygen enriched gas by adsorbing nitrogen to remove the 
nitrogen from the air; and 

an ultrasonic apparatus for measuring the 
concentration of the oxygen in the oxygen enriched gas 

3 0 and flow rate of the oxygen enriched gas, the ultrasonic 
apparatus comprising : 

a conduit for flowing an objective gas, the 
concentration of which is to be measured; 

a first ultrasonic transmission-reception device 

35 mounted to the inside of the conduit; 

a second ultrasonic transmission-reception device 
mounted to the inside of the conduit to face the first 
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ultrasonic transmission-reception device; 

the conduit includes a straight potion and 
perpendicular portions perpendicularly connected to the 
ends of the straight portion; 
5 the first and second ultrasonic transmission- 

reception devices are disposed in the perpendicular 
portions to face the ends of the straight portion; and 

the distance between the first and second ultrasonic 
transmission-reception devices and the respective ends of 
10 the straight portion of the conduit satisfying the 
following relation, 

0<D<f-r 2 /C 

where : 

D: the distance (m) between the first and second 
15 ultrasonic transmission-reception devices and the 
respective ends of the straight portion 
f: frequency of the ultrasonic waves in the sample gas 
(Hz) 

r: inner radius of the conduit (m) 
20 C: velocity of the ultrasonic waves (m/sec) 

Brief Description of the Drawings 

Figure 1 is a schematic diagram of an oxygen 
concentration apparatus according the invention; 
25 Figure 2 is a schematic diagram of an ultrasonic 

apparatus of the invention; 

Figure 3A is a waveform based on the received 
ultrasonic waves; 

Figure 3B is an enlarged illustration of a portion 
30 of the waveform shown in Figure 3A; 

Figure 4 is an illustration showing the ultrasonic 
waveform with the trigger signals and the zero-cross 
signals ; 

Figure 5 is a graph showing the relation between the 
35 velocity of the ultrasonic waves and the temperature; 

Figure 6 is an illustration showing forward and 
backward ultrasonic waveforms in case that the phases, at 
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which the trigger signals are generated, coincide with 
each other; 

Figure 7 is an illustration similar to that of 
Figure 6 in case that the phases do not coincide with 
5 each other; 

Figure 8 is an illustration similar to that of 
Figure 6 in case that the phases do not coincide with 
each other; 

Figure 9 is an explanatory illustration for 
10 explaining the way for obtaining the zero-cross time 
instant with an assumption that the sample gas is in 
stationary state; 

Figure 10 is an explanatory illustration for 
explaining the way for obtaining the ultrasonic reception 
15 point; 

Figure 11 is a section of the ultrasonic apparatus 
according to another embodiment of the invention; 

Figure 12 is an explanatory illustration for 
explaining the sound field formed at the front of an 
20 ultrasonic transducer; 

Figures 13 shows experimental results of ultrasonic 
waveforms which were obtained by an apparatus of Figure 
11; 

Figures 14 shows experimental results of ultrasonic 
25 waveforms which were obtained by an apparatus of Figure 
1 1 ; and 

Figures 15 shows experimental results of ultrasonic 
waveforms which were obtained by an apparatus of Figure 
11. 

30 

Best Mode for Carrying out the Invention 

A preferred embodiment of the present invention will 
be described below. In the embodiment described below, a 
sample gas is composed of a mixture of oxygen and 
35 nitrogen. However, the measurable sample gas is not 

limited to a sample gas of oxygen and nitrogen and the 
present invention can be supplied to a mixture including 
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another gas . 

Figure 1 shows a schematic diagram of an oxygen 
concentration system having an ultrasonic gas 
concentration and flow rate measuring apparatus according 
5 to a preferred embodiment of the present invention. 

The apparatus 100 includes an oxygen concentration 
apparatus 102 which produces an oxygen enriched gas by 
removing nitrogen from the air supplied by a compressor 
104 from the outside of the system through a filter 106. 

10 The oxygen enriched gas produced by the oxygen 

concentration apparatus 102 is supplied to an ultrasonic 
apparatus 200 of the present invention through a flow 
rate setting device 108, such as a pressure reduction 
valve. The produced oxygen enriched gas is then supplied 

15 to a user or a patient through a production filter 110. 

The oxygen concentration apparatus includes a 
plurality of columns (not shown) for accommodating 
nitrogen adsorbent such as a zeolite, a piping system 
(not shown) including conduits for directing the 

20 compressed air from the compressor 104. to each of the 
plurality of columns and for directing the produced 
oxygen enriched gas from the columns to the flow rate 
setting device 108 and valves (not shown) disposed in the 
piping system for selectively opening and closing the 

25 conduits so that the adsorbent contained in one of the 
columns adsorbs nitrogen to produce the oxygen enriched 
gas and the adsorbent contained in the other columns 
release the adsorbed nitrogen for the regeneration of the 
adsorbent. 

30 With reference to Figure 2, the ultrasonic apparatus 

200 of the present invention for measuring the 
concentration and the flow rate of a sample gas will be 
described below. 

The gas concentration and flow rate measuring 

35 apparatus 200 includes a conduit 202 for flowing a sample 
gas or the oxygen enriched gas produced by the oxygen 
concentration apparatus 102. The conduit 202 has a 
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straight portion 208 and perpendicular portions 204 and 
206 connected to the ends of the straight portion. The 
straight portion 208 comprises a conduit member having a 
circular section, the diameter of which does not changes 
5 along the longitudinal axis. A first ultrasonic 
transducer 218 , providing a first ultrasonic 
transmission-reception device, is fixedly provided at an 
end of the inside of the straight portion, and a second 
ultrasonic transducer 222, providing a second ultrasonic 

10 transmission-reception device, is fixedly mounted to the 
other end of the inside of the straight portion 208 to 
face the first ultrasonic transducer 218. In this 
embodiment, the distance between the first and second 
ultrasonic transducers 218 and 222 is referred to a 

15 propagation length L a . 

The perpendicular portion 204, disposed at the 
upstream side relative to the flow direction of the gas 
through the conduit 202, has an inlet port 204a. The 
oxygen concentration apparatus 102 is connected to the 

20 inlet port 204a as a sample gas source 212 through a 
supply conduit 210. 

The perpendicular portion 206, disposed at the 
downstream side relative to the flow direction of the gas 
through the conduit 202, has an outlet port 206a to which 

25 the production filter 110 is connected. 

A transmission-reception switch 224 is connected to 
the first and second ultrasonic transducers 218 and 2 22. 
The transmission-reception switch 224 independently 
switches the operation mode of the first and second 

30 ultrasonic transducers 218 and 222 between a transmission 
mode in which the first and second ultrasonic transducers 
218 and 222 transmit ultrasonic waves and a reception 
mode in which the first and second ultrasonic transducers 
218 and 222 receive the ultrasonic waves. The 

35 transmission-reception switch 224 is connected to a 
microcomputer 226 so that the switching operation of 
transmission-reception switch 224 is controlled by the 
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microcomputer 226. 

Temperature sensors 216 and 220, for measuring the 
temperature of the gas flowing through the conduit 202, 
are disposed preferably in the perpendicular portions 204 
5 and 206 so that they do not disturb the flow in the 

straight portion 2 08. The temperature sensors 216 and 
220 are connected to the microcomputer 226. In this 
connection, if the changes in the temperature of the 
sample gas is small, only one of the temperature sensors 

10 216 or 220 may be disposed. 

A driver 228 for driving the first and second 
ultrasonic transducers 218 and 222, a zero-cross 
detection circuit 230 for detecting zero-cross time 
instants of the signals from the first and second 

15 ultrasonic transducers 218 and 222, a display unit 234 
for indicating, for example, the operating condition of 
the device 200 and the measurement results and memory 232 
including a nonvolatile memory device or a disc device 
for storing the operation system for the microcomputer 

20 226 and various parameters are connected to the 

microcomputer 22 6 . 

The operation of the ultrasonic concentration and 
flow rate measuring apparatus 200 of the present 
embodiment will be described below. 

2 5 A sample gas, for example an oxygen-nitrogen gas 

mixture the mixture ratio of which is P: ( 1-P) ( OssPsl ) , is 
supplied to the conduit 202. At that time, the 
temperatures of the sample gas are measured by the 
temperature sensors 216 and 220 and the mean value 

3 0 thereof is stored in the memory 232 as a reference 

temperature T 0 (K). According to the embodiment, the 
working temperature range of the system 100 is preferably 
selected, for example 5-35 Celsius degrees. 

During the supply of the sample gas, pulses for 
35 generating the ultrasonic waves are transmitted to the 
driver 228 from the microcomputer 226. A pulse voltage 
is supplied to the first ultrasonic transducer 218 from 
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the driver 228 through the transmission-reception switch 
22 4. The first ultrasonic transducer 218 generates 
ultrasonic waves corresponding to the pulse voltage. The 
ultrasonic waves generated by the first ultrasonic 
5 transducer 218 propagate through the sample gas flowing 
through the straight portion 208 of the conduit 202 and 
are received by the second ultrasonic transducer 222. 
The second ultrasonic transducer 222 generates an 
electric signal corresponding to the received ultrasonic 

10 waves to the microcomputer 22 6 through the transmission- 
reception switch 224 and the zero-cross detection circuit 
230. The microcomputer 226 calculates the forward 
propagation time t al (sec) on the basis of the time when 
the transmitted pulses are generated to the driver 228 

15 and the time when the electric signal is received from 
the second ultrasonic transducer 222. 

Then, the transmission-reception switch 224 switches 
the operation mode of the first ultrasonic transducer 218 
from the transmission mode to the reception mode right 

20 after the electric signal from the second ultrasonic 
transducer 222 is received and also switches the 
operation mode of the second ultrasonic transducer 222 
from the reception mode to the transmission mode. 
Thereafter, pulses for generating the ultrasonic waves 

25 are transmitted to the driver 228 from the microcomputer 
226. A pulse voltage is supplied to the second 
ultrasonic transducer 222 from the driver 228 through the 
transmission-reception switch 224. The second ultrasonic 
transducer 22 2 generates ultrasonic waves corresponding 

30 to the pulse voltage. The ultrasonic waves are received 
by the first ultrasonic transducer 218. The first 
ultrasonic transducer 218 generates an electric signal 
corresponding to the received ultrasonic waves to the 
microcomputer 226 through the transmission-reception 

35 switch 224 and the zero-cross detection circuit 230. The 
microcomputer 22 6 calculates the backward propagation 
time t a2 (sec) on the basis of the time when the 
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transmitted pulses are generated to the driver 228 and 
the time when the electric signal is received from the 
first ultrasonic transducer 218, 

By obtaining the mean value of t 8l and t s2 , the 
5 affection of the flow of the sample gas in the conduit 
202 can be removed . The ultrasonic propagation time t s 
in the stationary sample gas is defined by following 
equation (4 ) . 

t 8 =(t sl + t 82 )/2 ...(4) 
10 The microcomputer 226 then calculates the ultrasonic 

propagation velocity C s (m/sec) through the stationary 
sample gas by flowing equation (5). 

C s =L s /t s .(5) 
The concentration of oxygen P s is obtained by 
15 following equation (6) on the basis of equations (1) and 
(2). 

P S =(KRT S /C S 2 -M H2 )/(M 02 -M N2 ) ...(6) 
Further, the concentration of oxygen in the sample 
can be obtained as a ratio of the ultrasonic propagation 

20 velocity in the sample gas and the ultrasonic propagation 
velocities in 100% of oxygen gas and 100% of nitrogen 
gas. That is, the ultrasonic propagation velocity 
C 02 (m/sec) at temperature T S (K) through 100% of oxygen gas 
and the ultrasonic propagation velocity C N2 (m/sec) at 

25 temperature T S (K) through 100% of nitrogen gas can be 

easily obtained by using equation (1). Thus, P s can be 
calculated by following equation (7) with the ultrasonic 
propagation velocity C s (m/sec) through the sample gas. 
P s = ( 1 /C s 2 - 1 /C N2 2 ) / ( 1 /C G2 2 - 1 /C N2 2 ) . . . ( 7 ) 

30 Such calculations are conducted by the microcomputer 

126, and the results are indicated by the display unit 
134. 

Next, the explanation will be directed to a method 
of obtaining t 8l and t a2 . In this connection, the moment 
35 when the first or second ultrasonic transducer 218 or 222 
transmits the ultrasonic waves is referred to as an 



WO 2004/094960 



- 16 - 



PCT/JP2004/005590 



emission time and the moment when the first or second 
ultrasonic transducer 218 or 222 receives the ultrasonic 
waves is referred to as an ultrasonic reception point in 
the present application. 
5 Figure 3A shows a typical ultrasonic waveform 

received by the microcomputer 22 6 and Figure 3B is an 
enlargement of a portion of the waveform shown by circle 
3B. As shown by Figures 3A and 3B, the waveform includes 
a various noise components which make difficult to detect 

10 the ultrasonic reception point of the ultrasonic waves 

propagating through the sample gas. Therefore, according 
to the present invention, the ultrasonic reception point 
is estimated on basis of the zero-cross time instant of 
the waveform which is detected after the amplitude of the 

15 waveform sufficiently increases to an extent. For this 
purpose, the zero-cross detection circuit 230 includes a 
zero-cross comparator and a trigger comparator. 

Referring to Figure 4 , the trigger comparator 
outputs a trigger signals S ti to the microcomputer 226 

20 when the waveform upwardly passes over a predetermined 
level. The zero-cross comparator outputs a zero-cross 
signal Z ci to the microcomputer 226 when the waveform 
upwardly passes over the zero level. The microcomputer 
226 determines each of the zero-cross signals Z ci as a 

25 zero-cross time instant after the first trigger signal 

S tl is received by the microcomputer 226. Preferably, 
the microcomputer 22 6 determines first three zero-cross 
signals as first to third zero-cross time instants Z cl , 
Z c2 and Z c3 . 

30 The interval between each of the zero-cross time 

instants theoretically corresponds to the cycle of the 
ultrasonic waves. Therefore, the ultrasonic reception 
point can be estimated by tracing back from the fist 
zero-cross time instant. Z cl along the time axis by an 

35 integral multiple of the cycle of the ultrasonic waves 
and therefore, the propagation time can be estimated by 
subtracting the emission time and an integral multiple of 
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the cycle of the ultrasonic waves from the ultrasonic 
reception point. 

As described above , the velocity of ultrasonic waves 
C(m/sec) propagating through a stationary gas is 
5 presented by equation (1). For example, the velocity of 
ultrasonic waves through a pure nitrogen gas at 2 0 
Celsius degrees is 349.1 m/sec, and the velocity of 
ultrasonic waves through a pure oxygen gas at 20 Celsius 
degrees is 326.6 m/sec. Therefore, at 20 Celsius 

10 degrees, the velocity of ultrasonic waves through an 
oxygen-nitrogen gas mixture falls within the rage of 
326.6 to 49.1 m/sec. Figure 5 is a graph showing the 
relation between the velocity of ultrasonic waves and the 
gas temperature, in which the upper and lower limits of 

15 the ultrasonic velocity through an oxygen-nitrogen gas 

mixture are indicated by C^fT) and C^CT). The possible 
propagation time rage is I^/C^fT) to L,/C ain (T). 
Therefore, if the propagation length L s is selected to 
satisfy following relation (8), only one integer can be 

2 0 selected which allows the ultrasonic reception point to 
fall within the possible propagation time range. 

(Ls/CWT) - L./C^fT)) < 1/f ...(8) 

where : 

f: frequency of the ultrasonic waves in the sample gas 
25 The gas temperature T which gives the maximum value 

of (Lg/C^fT) - L e /C inin (T)) is the lower limit of the 
working temperature. If the working temperature is 5 
Celsius degrees, and the frequency of the ultrasonic 
waves is 40 KHz , the propagation length L e which 
30 satisfies relation (8) is calculated as follows. 

L s < 12.3 cm ...(9) 

According to the embodiment, L g =0.1 m is employed as 
an example. 

In order to obtain the ultrasonic propagation time 
35 t s , the forward and backward propagation times t sl and t o2 
is previously measured. With reference to Figure 6, the 
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trigger signals are generated when the second waves in 
both the forward and backward waveforms passes over the 
trigger level . In this case, the trigger signals are 
generated at the same timing or phase relative to the 
5 waveforms and the difference in the zero-cross time 

instants between the forward and backward waves, A = Z cBi 
- Z cFi , is substantially equal to the difference t d in 
propagation times t sl and t s2 between the forward and 
backward waves (Z cFi : the zero-cross time instants of the 

10 forward waveform, Z cBi : the zero-cross time instants of 
the backward waveform, i=l, 2, 3... (wave number)). 

However, the trigger signals S ti are often generated 
at different phases of the waveforms between the forward 
and backward waves, even if the same trigger level is 

15 used. Referring to Figure 7, for the forward waves, the 
trigger signal is generated when the third wave passes 
over the trigger level and for the backward waves, the 
trigger signal is generated when the second wave passes 
over the trigger level. Therefore, the trigger signal 

20 for the backward waves is generated at one cycle earlier 
than the trigger signal for the forward waves. In this 
case, the difference in the zero-cross time instants 
between the forward and backward waves, A = Z cBi - Z cFi , 
gives a negative value. If the sample gas flows through 

25 the conduit 102, A = Z cBi - Z cFi must not become negative. 

Therefore, if A = Z cBi - Z cFi gives a negative value, it is 
apparent that the trigger signal for the backward waves 
is generated earlier than the trigger signal for the 
forward waves . 

30 On the other hand, referring to Figure 8, for the 

forward waves, the trigger signal is generated when the 
second wave passes over the trigger level and for the 
backward waves, the trigger signal is generated when the 
third wave passes over the trigger level. In this case, 

35 the difference in the zero-cross time instants between 
the forward and backward waves, A = Z cBi - Z cFi , exceeds 
one cycle of the ultrasonic waves, which indicates that 
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the trigger signal for the forward waves is generated 
earlier than the trigger signal for the backward waves • 
According to the embodiment of the invention, the 
conduit 102 is designed so that the propagation time 
5 difference t d between the forward and backward waves 

always falls within one cycle of the ultrasonic waves. 
This feature allows the microcomputer 226 to distinguish 
the cases shown in Figures 7 and 8 from each other and to 
calculate the propagation time difference t d . That is, 
10 if A=Z cBi -Z cFi is negative, the case is as shown in Figure 
7 and if A = Z^-Z^ exceeds one cycle of the ultrasonic 
waves the case is as shown in Figure 8 . 

Thus, the configuration of conduit 102 which has the 
above feature will be described below. 
15 The possible range of the flow velocity V (m/sec) of 

the sample gas is presented by following inequality (10). 

0<;Vs;Q/(600007T/r 2 ) . . (10) 

where: 

Q: flow rate of the sample gas (litter/min) 

20 r: inner radius of the conduit (cm) 

As described above, the velocity of ultrasonic waves 
propagated in the forward direction relative to the 
sample gas flow is C^C+V, and the velocity of ultrasonic 
waves propagated in the backward direction relative to 

25 the sample gas flow is C 2 =C-V. 
Where : 

C: the velocity of ultrasonic waves propagating through a 
stationary sample gas (m/sec) 

C x : the velocity of the ultrasonic waves propagated in 
3 0 the forward direction relative to the sample gas flow 
(m/sec) 

C 2 : the velocity of ultrasonic waves propagated in the 
backward direction relative to the sample gas flow 
(m/sec) 

35 V: the flow velocity (m/sec) 

The propagation time difference t d is calculated by 
the following equation. 
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t d - L 8 /C 2 -L 8 /C 1 

- L./(C-V)-L./(C+V) ...(11) 
Therefore, if the inner radius of the conduit 102 
satisfies following relation (12), the propagation time 
5 difference t d becomes smaller than the cycle of the 
ultrasonic waves. 

L 8 /(C-Q/(60000jt/r 2 ) )- L 8 / (C+Q/ ( 60000ji/r 2 ) )<l/f ...(12) 
The left term of inequality (12) is maximized when 
the velocity of the ultrasonic waves through the conduit 
10 102 is minimum (0=0^(5 Celsius degrees ) =318 . 1 m/sec). 

Therefore, for example, if the frequency of the 
ultrasonic waves through the conduit 102 is 40 (KHz), the 
flow rate Q=10 ( litter /min) and the length of the conduit 
102 is 10 (cm), then the inner radius r (mm) of the 
15 conduit 102 is r>2.05 (mm)* According to the embodiment, 
r=2.5 (mm) is selected as an example* 

Next, the method for measuring the concentration and 
flow rate of a sample gas will be described in detail 
below. 

20 First, in case shown in Figure 6, the propagation 

time difference t d between the forward and backward waves 
is obtained by A = Z cBi - Z cFi , because, as described 
above, the propagation time difference t d is 
substantially equal to the difference A = Z cBi - Z cFi . In 

25 case shown in Figure 7, the propagation time difference 
t d is obtained by B = Z cBi ^ - Z cFi . Further, in case shown 
in Figure 8, the propagation time difference t d is 
. obtained by B = Z cBi - Z cFi+1 . Preferably, a plurality of 
the values A or B are obtained for arithmetic average. 

30 Next, the velocity of the ultrasonic waves through 

the sample gas is estimated with the assumption that the 
sample gas is in stationary state. For this purpose, the 
phase difference in the outputs of the trigger signals is 
previously determined on the basis of the value of A. If 

35 there is no phase difference, as shown in Figure 6, mean 
value Z c ave of first zero-cross time instants of the 
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forward and backward waveforms is calculated by the 
following equation* 

Z c .ave=(Zcri + Z cB1 )/2 ...(13) 

In case shown in Figure 7, mean value Z cave of first 
5 zero-cross time instants of the forward and backward 
waveforms is calculated by the following equation. 

Zc.ave=(ZcFl + Z cB2 )/2 ...(14) 

In case shown in Figure 8, mean value Z cave of first 
zero-cross time instants of the forward and backward 
10 waveforms is calculated by the following equation, 

Z c _ave=(Z cF2 + Z cB1 )/2 ...(15) 

The mean value Z cave can be considered as the first 
zero-cross time instant which is obtained with the 
assumption of the ultrasonic waves through the stationary 

15 sample gas. Z c avo is referred to a reference zero-cross 
time instant in the present application. 

As described above, the length of the conduit 102 is 
designed so that only one integer is selected allowing 
the ultrasonic reception point to fall in the possible 

20 range of the propagation time (Figure 9). Therefore, the 
ultrasonic propagation time t 8 is estimated by tracing 
back from the fist zero-cross time instant Z c ave along the 
time axis by an integral multiple of the cycle of the 
ultrasonic waves until the ultrasonic reception point 

25 falls in the possible range. 

The velocity of the ultrasonic waves Cs through the 
stationary sample gas is estimated by following equation 
(16). 

C 8 = L./t. ...(16) 
30 The concentration of oxygen P s is obtained by 

equation (6) or (7) with the calculated C e . 

The forward and backward propagation times t sl and 
t a2 through the sample gas flowing through the conduit 
102 are estimated by following equations (17) and (18). 
35 t eX =t s -t d /2 ...(17) 
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t 8l =t s +t d /2 ...(18) 
The forward and backward velocities C x and C 2 of the 
ultrasonic waves through the sample gas flowing through 
the conduit 102 are estimated by following equations (19) 
5 and (20). 

Cx=L s /t 8l ...(19) 
C 2 =L 8 /t s2 ...(20) 
Then, the flow velocity V of the sample gas through 
the conduit 102 is obtained by equations (3), (19) and 
10 (20). Further, the flow rate Q of the sample gas is 
calculated by following equation (21). 

Q=6000jrr 2 V ...(21) . 
Next, with reference to Figures 11 to 15, a 
preferred embodiment of the ultrasonic concentration and 
15 flow rate measurement apparatus will be described below, 
the ultrasonic concentration and flow rate 
measurement apparatus 10 includes a conduit 27 which 
provides the conduit 102 of the embodiment of Figure 2. 
Housings 25 and 26, for enclosing first and second 
20 ultrasonic transducers 20 and 21, secured to the ends of 
the conduit 2 7 by welded portions 41 and 42. The 
housings 25 and 26 include ports 28 and 29 extending 
perpendicular to the conduit 27 to provide the inlet and 
outlet portions 204a and 206a of the embodiment of Figure 
25 2. The conduit 27 and the housings 25 and 26 are 

preferably made of the same metallic material such as an 
aluminum alloy. 

The conduit 27 and the housings 25 and 26 are 
secured at one point to a substrate 30 or a housing of 
30 the oxygen concentration apparatus by a screw 45. This 
configuration allows the longitudinal deformation of the 
conduit 27 freely from external force which may be 
generated when the conduit 2 7 is thermally deformed. 

Covers 23 and 24 are attached to the housings 25 and 
35 26 to close the end openings of the housings by screws 43 
and 44 with O-rings 39 and 40 being clamped between the 
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housings 25 and 26 and the covers 23 and 24. The first 
and second ultrasonic transducers 20 and 21 are attached 
to the inner surfaces of the covers 23 and 24. The first 
and second ultrasonic transducers 20 and 21 generate 40 
5 KHz of ultrasonic waves. 

Further, temperature sensors 37 and 38 for detecting 
the gas temperature are attached to the inner surfaces of 
the covers 23 and 24. The first and second ultrasonic 
transducers 20 and 21 and the temperature sensors 37 and 

10 38 are coupled to the microcomputer 226 through 

connectors 31 and 34 attached to the outer surfaces of 
the covers 23 and 24 , cables 33 and 36 and connectors 32 
and 35 mounted on the substrate 30. 

Distance D between the end faces of the first and 

15 second ultrasonic transducers 20 and 21 and the 

respective ends of the conduit 27 is a significant design 
matter. Generally, the sound field formed by the 
ultrasonic waves from an ultrasonic transducer includes 
near and far sound fields, as shown in Figure 12. The 

20 ultrasonic waves propagate linearly through the near 
sound field and on the other hand, in the far sound 
field, spread in the form of spherical waves. Therefore, 
if the ends of the conduit 27 are out of the near sound 
field, the ultrasonic energy transmitted in the conduit 

25 2 7 is reduced compared with conduit having the ends 
disposed in the near sound field and therefore, the 
sound/noise ratio of the signal from the transducers is 
reduced . 

It is known that the boundary between the near and 
30 far sound fields is presented at a point Z 0 the distance 
D of which from the end face of an ultrasonic transducer 
along the center line of the transducer is defined by 
following equation (22). 

D=f-r 2 /C ...(22) 

35 Where: 

f : frequency of the ultrasonic waves in the sample gas 
(Hz) 
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r: inner radius of the conduit (m) 

C: velocity of the ultrasonic waves (m/sec) 

As described above, the velocity C through a sample 
gas is defined by equation (1). Therefore, the higher 
5 the gas temperature and the smaller the molecular weight, 
the higher the velocity C becomes. According the 
embodiment, the condition which maximizes Z 0 is that, for 
example, the sample gas is the air at 35 Celsius degrees 
and then, Z 0 is about 1.4 mm. 

10 Figures 13-15 show experimental results of 

ultrasonic waveforms which were obtained by an apparatus 
of Figure 11 with the distance d of 0.3 mm, 1.0 mm and 
1.8 mm. The experimental results suggest that the 
ultrasonic energy received by the ultrasonic transducer 

15 is significantly reduced when the distance d is 1.8 mm 
compared with the cases of the distance d of 0.3 mm and 
1.0 mm. 



